Passage Times for Polymer Translocation Pulled through a Narrow Pore 
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We study the passage times of a translocating polymer of length N in three dimensions, while it 
is pulled through a narrow pore with a constant force F applied to one end of the polymer. At 
small to moderate forces, satisfying the condition FN" /kgT < 1, where v ~ 0.588 is the Flory 
exponent for the polymer, we find that tjv, the mean time the polymer takes to leave the pore, 
scales as jV 2 +" independent of F, in agreement with our earlier result for F = 0. At strong forces, 
i.e., for FN" I kgT 3> 1, the behaviour of the passage time crosses over to rjv ~ N 2 /F. We show 
here that these behaviours stem from the polymer dynamics at the immediate vicinity of the pore 
— in particular, the memory effects in the polymer chain tension imbalance across the pore. 
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[. INTRODUCTION 

Molecular transport through cell membranes is an es- 
sential mechanism in living organisms. Often, the mole- 
cules are too long, and the pores in the membranes too 
narrow, to allow the molecules to pass through as a single 
unit. In such circumstances, the molecules have to de- 
form themselves in order to squeeze — i.e., translocate — 
themselves through the pores. DNA, RNA and proteins 
are such naturally occuring long molecules U; UG| 3 13) 
in a variety of biologicalprocesses. Translocation is also 
used in gene therapy f^; 0) , in delivery of drug molecules 
to their activation sites (@) , and as a potentially cheaper 
alternative for single-molecule DNA or RNA sequencing 
|1; [IB). Consequently, the study of translocation is an 
active field of research: as a cornerstone of many biolog- 
ical processes, and also due to its relevance for practical 
applications. 

Translocation in living organisms is a complex process. 
Take for instance the case of gene expression: most pro- 
teins are synthesized within the cytoplasm. Their subse- 
quent accurate and swift delivery to target sites, requir- 
ing energy, is a crucial step in gene expression. In differ- 
ent situations the energy is provided by chaperon mole- 
cules 

GH), pHgradient (0) or molecular motors across 
membranes lfl3l ) . These delivery mechanisms can be fur- 
ther complicated by membrane fluctuations and some- 
times by gates that control the accessibility of the pores 
. In view of such complexity, translocation as a bio- 
logical or biophysical process in living organisms has been 
scrutinized in a variety of in vivo experimental situations 
[see e.g. Ref. (fl5h and the references therein]. 

The experimental developments have been followed by 
a number of mean-field type theoretical studies on poly- 
mer translocation 1|2(J). 

More recently, translocation has found itself at the 
forefront of single-molecule-detection experiments |f| 
[l7h . as new developments in design and fabrication of 
nanometer-sized pores and etching methods may lead to 



cheaper and faster technology for the analysis and detec- 
tion of single macromolecules. The underlying principle 
for these experiments is that of a Coulter counter: mole- 
cules suspended in an electrolyte solution pass through a 
narrow pore in a membrane. The electrical impedance 
of the pore increases with the entrance of a molecule 
as it displaces its own volume of the electrolyte solu- 
tion. By applying a voltage over the pore, the passing 
molecules are detected as current dips. For nanometer- 
sized pores (slightly larger than the molecule's cross- 
section) the magnitude and the duration of these dips 
have proved to be effective in determining the size and 
length of the molecules. In the case of DNA sequencing 
at nucleotide level, usage of protein pores (modified a- 
haemolysin, mitochondrial ion channel, nucleic acid bind- 
ing/channel protein etc.), and etching specific DNA se- 
quences inside the pores @; [T3) have opened up promis- 
ing new avenues of fast, simple and cheap technology for 
single macromolecule detection, analysis and characteri- 
zation [see Ref. lfl9l ) for a recent development]. 

The subject of this paper is a translocating polymer 
threaded through a narrow pore in an immobile mem- 
brane, where a bead is attached to one end of the poly- 
mer, and the bead is pulled by an optical tweezer with a 
constant force. Such a setup can be used to spread apart 
a partially unzipped dsDNA molecule — of which one 
strand is threaded through the pore — a process that 
can quantify the forces involved in basepair unzipping 
kinetics l|2lh. In theoretical literature, this problem has 
been considered in recent times: for polymer length N 
and applied force F, in Ref. H3), in the absence of hy- 
drodynamical interactions, a lower bound oc N 2 for small 
forces (FN U /ksT < 1) has been argued for the polymer's 
mean unthreading time tm, the average time it takes for 
the polymer to leave the pore. The lower bound holds 
in the limit of unimpeded polymer movement, i.e., for an 
infinite pore, or equivalently, in the absence of the mem- 
brane. Simulation data (in two-dimensions) presented in 
Ref. ((22h indicated that the lower bound may very well 
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be valid in the limit of narrow pores as well. The same 
problem, also in the absence of hydrodynamical interac- 
tions, has been numerically studied in two-dimensions in 
Ref. (23] ). It reported that for narrow pores tn ~ iV 2 
with the velocity of translocation v{t) ~ TV -1 for moder- 
ate and strong forces. As the force-dependence of tn is 
concerned, Ref. HI) reported numerical results that in 
the absence of the membrane tn ~ F~ 2+ ^ for moderate 
forces, and tn ~ F^ 1 for strong forces, while for narrow 
pores tn ~ F^ 1 for moderate to strong forces. 

The purpose of this paper is to revisit the problem 
of translocation of a polymer pulled through a narrow 
pore in the absence of hydrodynamical interactions, in 
order to provide a deeper theoretical understanding of 
the polymer dynamics under these conditions, as well as 
of the scaling behaviour of the unthreading time of the 
polymer. In support of our theory, we perform high pre- 
cision computer simulations, using a three-dimensional 
self-avoiding lattice polymer model that we have used 
before to study polymer translocation H3; HE; H3) and 
several other situations |28h . Our conventions to study 
this problem, all throughout this paper, is the following. 
We place the membrane at z — 0, and thread the poly- 
mer of total length 2N halfway through the pore such 
that both the right (z > 0) and the left (z < 0) of the 
membrane have equal number of monomers N. We fix 
the middle monomer (monomer number N) at the pore, 
apply a force F on the free end on the right and let both 
left and right segments of the polymer come to equilib- 
rium. At t — we release the middle monomer and let 
translocation commence. The mean time tn that the 
polymer remains within the pore is defined as the mean 
unthreading time for polymer length N under the force 
F. Additionally, we use fcsT = 1, although fc^T is ex- 
plicitly mentioned at several places in the paper. 

Our main results in this paper are as follows. At small 
to moderate forces, satisfying the condition FN" < 1, 
where v w 0.588 is the Flory exponent for the poly- 
mer, we find that tn is independent of F . In agreement 
with our earlier result for unbiased polymer transloca- 
tion; i.e., for F — (24), tn scales with polymer length 
as t n ~ N 2+v . At strong forces, i.e., for FN" > 1, we 
find tn ~ N 2 /F . While these results agree with the ex- 
isting ones |22l :[23h in broad terms, we show that v(t), the 
velocity of translocation is not constant in time. In fact, 
for strong forces, we show that the velocity of transloca- 
tion v(t) behaves as i -1 / 2 , while for small to moderate 
forces the behaviour of v(t) is more complicated. The 
physical picture provided in Refs. |22l; [23). wherein the 
scaling arguments for the unthreading time involved a 
constant velocity of translocation (albeit an average one, 
in light of this work) is incomplete |29t ) . Using theoreti- 
cal analysis supported by high-precision simulation data, 
we show that these behaviours stem from the dynamics 
of the polymer segments at the immediate vicinity of the 
pore — in particular, the memory effects in the polymer 
chain tension imbalance across the pore. The theoretical 
analysis presented here is based on that of Ref. ((24T 1 , and 



therefore provides a direct confirmation of the robustness 
of the theoretical method presented in Ref. |24h . 

This paper is organized in the following manner. In 
Sec. [n] we discuss a method to measure component of 
the polymer chain tension which is perpendicular to the 
membrane. In Sec. Mil we analyze the memory effects 
in 4>(t), the imbalance of this component of the polymer 
chain tension. In Sec. IIVI we discuss the consequence of 
these memory effects on the translocation velocity v(t), 
and obtain the relation between the mean unthreading 
time tn and the polymer length N. We finally end this 
paper with a discussion in Sec. LY1 

II. CHAIN TENSION PERPENDICULAR TO THE 
MEMBRANE 

A translocating polymer should be thought of as two 
segments of polymers threaded at the pore, while the 
segments are able to exchange monomers between them 
through the pore. In Ref. (|24T ) we developed a theoreti- 
cal method to relate the dynamics of translocation to the 
imbalance of chain tension between these two segments 
across the pore. The key idea behind this method is that 
the exchange of monomers across the pore responds to 
<fr(t), this imbalance of chain tension; in its turn, (p(t) ad- 
justs to v(t), the transport velocity of monomers across 
the pore. Here, v(t) = s(t) is the rate of exchange of 
monomers from one side to the other, where s(t) is the 
total number of monomers transferred from one side of 
the pore to the other in time [0,t], In fact, we noted that 
s(t) and 4>{t) are conjugate variables in the thermody- 
namic sense, with <f)(t) playing the role of the chemical 
potential difference across the pore. 

By definition, <f>(t) = $B.{t) ~ where and 

$l(£) are respectively the chain tension (or the chemi- 
cal potential) on the right and the left side of the pore. 
Consider a separate problem, where we tether one end 
of a polymer to a fixed membrane, yet the number of 
monomers are allowed to spontaneously enter or leave 
the tethered end, then we have 

Wt( > +) 

yfc - _j = eMHt)/k B T} , (1) 

where Wt{ ► +) [resp. Wt(-\ > — )] is the rate that a 

monomer enters (resp. leaves) the polymer chain through 
the tethered end at time t, Note that tethering the poly- 
mer while allowing monomers to enter or leave the poly- 
mer at the tethered end is precisely the case that trans- 
location represents. 

Returning to our problem of a translocating polymer 
under a pulling force F, note that at t = 0, when the 
left and the right segments are equilibrated with F = 
and F ^ respectively, it is easy to use Eq. JT]) to 
measure the chain tension for both segments at the pore 
[$(i = 0) in our notation], since under these conditions, 
we also have the relation that 

P-W t =o(- -> +) = P+Wt=o(+ -> -) , (2) 
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where P_ (resp. P+) is the probability that the (left or 
the right) polymer segment has one monomer less (resp. 
one extra monomer). Equations HJ and ^ together 
yield us 

$(t = 0) = fc B Tln|±. (3) 

Note that even for F = 0, there is nonzero chain ten- 
sion, due to the presence of the membrane. A polymer's 
free energy close to a membrane is higher than its free 
energy in bulk. In other words, the membrane repels 
the polymer, and as a result, for a polymer with one 
end tethered to a membrane, the monomers close to the 
membrane are more stretched than they would be in the 
bulk. 

The chain tension as obtained from Eq. ([3]) is linearly 
related to the z-coordinate of the centre-of-mass of the 
first few monomers along the polymer's backbone, at the 
immediate vicinity of the pore, at least for the relatively 
modest forces used in our simulations. This is shown in 
Fig. [TJ where for a tethered polymer of length N = 100, 
the average distance (Z^'(t — 0)} of the centre-of-mass 
of the first 4 monomers along the polymer's backbone, 
counting from the tethered end of a polymer, is plotted 
versus the chain tension <&, while its free end is pulled 
with various force strengths F. Within the error bars, all 
the points in Fig. [TJfall on a straight line, implying that 
<I> is very well-proxied by (Z^). Note in Fig. [TJthat the 
black line does not pass through the origin, which shows 
that $f=o 7^ 0, as we argued above. Since measurements 
of the chain tension via Eq. ([3]) are much more noisy than 
measurements of (Z^), we will use the latter quantity 
as a measure for the chain tension. 




1.3125 1.4375 1.5625 1.6875 



< T (t = 0)> 

FIG. 1 {Z (4) (t = 0)) vs. $(t = 0) demonstrating the lin- 
ear relationship between the two, for N = 100 and F = 
0.0,0.3,0.5,0.7,0.9 and 1.0 respectively. The angular brack- 
ets for {Z w (t = 0)) indicates an average over 12, 800, 000 
polymer realizations. The data for for &(t — 0) are obtained 
2, 400 polymer realizations. The red line corresponds to the 
linear best-fit. Inset: (Z^{t = 0)) as a function of F. 



III. MEMORY EFFECTS IN THE z-COMPONENT OF THE 
CHAIN TENSION 

In the case of unbiased polymer translocation, we 
have witnessed in Ref. (|24r i that the memory effects of 
the polymer gives rise to anomalous dynamics of trans- 
location. We argued j2~4T l that the imbalance of the 
chain tension <f>(t) across the pore and the number of 
monomers s(t) that have crossed from one side of the 
membrane to the other in time [0, t] are conjugate vari- 
ables in the thermodynamic sense. Additionally, <j)(t) is 
related to the translocation velocity v(t) by the relation 
4>{t) = 4>t—o — Jq dt'n(t — t')v(t') via the memory kernel 
^(t), which can be thought of as the 'impedance' of the 
system. On average, there will not be an imbalance in 
chain tension if no force is applied, but there will be fluc- 
tuations in chain tension. When the polymer is pulled by 
a force F to the right, the symmetry between the poly- 
mer segments on two sides of the membrane (viz., the 
polymer segments on the right are more stretched than 
those on the left of the membrane), are destroyed. As a 
result, on average <f>(t), 4>t=o and v(t) are non-zero, and 
from now on, we understand these three quantities as an 
average over all the unthreading polymers. Additionally, 
for F ^ the memory effects continue to be present, and 
the memory kernels ^h{t) and (J-n(t) for the polymer seg- 
ments on the left and the right sides of the membrane are 
different. In Ref. |24~1 ) we determined /ii(i) = fJ.F=o(t) 
by tethering a polymer of length N — 10 on a fixed mem- 
brane, where we injected p monomers at the tethered 
end at time t = 0, i.e., v(t) = pS(t) with p = 10 (bring- 
ing the final polymer length to N), and proxying <f>(t) by 
the average distance of the centre-of-mass of the first 5 
monomers {Z^(t)) from the membrane. We found 

~ exp(-t/r Rouse ) , (4) 

where TR OUSe ~ N 1+2u is the Rouse time, the longest 
relaxation time-scale of a polymer of length N. 

Following the same line as presented in Ref. ((241 ) , here 
we compute HB,{t), the memory effect of a polymer of 
length N with one end tethered to a membrane, and the 
other end pulled by a force F. The first step to do this 
is to obtain the relaxation time for a polymer of length 
N under these conditions. For F = the result for the 
relaxation time ~ TR 0use is well-known [and has been con- 
firmed in an earlier study of ours (27,)], but with F ^ 0, to 
the best of our knowledge, the corresponding analytical 
result does not exist. We therefore resort to simulations: 
we denote the vector distance of the free end of the poly- 
mer w.r.t. the tethered end at time t by e(t), and define 
the correlation coefficient for the end-to-end vector as 

c(t) = <e(t) ■ e(0)) - (e(t)) ■ (e(0)) 

V<e 2 W-(eW) 2 )(e 2 (0)-(e(0)) 2 )- 

The angular brackets in Eq. © denote simple ensemble 
averaging for F ^ 0. We first obtain the time correlation 
coefficients c(t) for 256 independent polymers, and c(t) 
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FIG. 2 (a) c(t) ~ exp ( — t j tf ) for strong forces, with tf ~ A'' 2 ; the data shown correspond to f 1 = 1.0; data obtained using 
256 polymers for each value of N. Inset: the same data are shown in semi-log plot to show that the decay of c(t) in time is 
exponential at long times, (b) Behaviour of fJ.R.(t), proxied by (Z' 4 '(t)), for N = 100 and four different values of F: F — 0.0 
(maroon), F = 0.3 (red), F = 0.5 (green), and F — 1.0 (blue); the solid black line corresponds to a slope £ -1//2 ; data obtained 
using 12, 800, 000 polymers for each value of F. See text for more details. 



is a further arithmatic mean of the corresponding 256 
different time correlation coefficients. At strong forces, 
when we scale the units of time by factors of N 2 (for self- 
avoiding polymers!), the c(t) vs. t curves collapse on top 
of each other. This is shown in Fig. [2ja) for F — 1.0 and 
N = 100, ...,350. 

What happens at small to moderate forces to the re- 
laxation time is not entirely clear to us. We do not expect 
the relaxation time-scale to change continuously with F. 
Thus, given the two limits r Rouse ~ N 1+2y for F = 
and tf ~ N 2 for strong forces, we believe that at small 
to moderate forces the relaxation time becomes a linear 
combination of TR 0use ~ 

N 1+2u and tf ~ N 2 , with the 
coefficients of these two times varying with the magni- 
tude of F. 

While Fig. [5Ja) provides the answer to the relaxation 
of the entire polymer for strong forces, the second step 
to identify a for nii{t) ~ t~ a exp(— t/rp) for some a for 
strong forces is to analyze the relaxation of the polymer 
segments at the immediate vicinity of the tethered point. 
The value of a depends on the relaxation properties fol- 
lowing the event of injecting, say, p extra monomers at 
the tether end, just like extra monomers add to (or get 
taken out of) the right segment of the polymer dur- 
ing translocation. Given the exp(— £/tf) behaviour of 
Fig. Eta), we anticipate that by time t after the extra 
monomers are injected at the tethered point, the extra 
monomers will come to a steady state across the inner 
part of the polymer up to n t ~ t 1 ! 2 monomers from the 
tethered point, but not significantly further. This inter- 
nal section of n t + p monomers in steady state extends 
only to r(n t ) from the membrane, because the larger scale 
conformation has yet to adjust, and consequently there 



is a compressive force / on these nt monomers. 

For F = 0, r{n t ) ~ n v t is the only length scale for 
the equilibrated section of the chain, which leads to 
f F =o ~ k B T5r/r 2 Hi), but for F ^ this does not hold. 
For small to moderate forces, i.e., for FN V /ksT < 1, the 
polymer conformation is given by a sequence of blobs of 
size £, given by the relation F£ = k B T and for strong 
forces, i.e., for FN u /k B T S> — ► a, where a is the 
size of a single monomer {H; I25T I. Thus, for F ^ 0, 
the shape of the polymer resembles that of a cylinder, 
implying /f^o ~ kB,T8r/{r£). The independence of £ 
on n t implies that r(n t ) ~ n t , which allows us to write 
/f^o ~ ksTSnt/ '(ntC) ~ t" 1 ^ 2 . This force is transmitted 
to the membrane, through a combination of decreased 
tension at the tether and increased incidence of other 
membrane contacts. The fraction borne by reducing ten- 
sion leads us to what is, strictly speaking, an inequality: 
a > 1/2. However, it seems unlikely that the adjust- 
ment at the membrane should be disproportionately dis- 
tributed between the two nearly balancing effects of poly- 
mer chain tension and monomeric repulsion, leading to 
the expectation that the inequality becomes an equality. 

Theoretically however, we cannot rule out the larger 
values for a, but our numerical results in Fig. E^b), 
where we have used (Z^\t)) to proxy &r(£) and p = 5, 
for strong forces favour the smallest theoretical value, 
namely a = 1/2. The power law decay preceding the ex- 
ponential ones in Fig. [2jb) change from t~ x + 2 » at F = 
tof- 1/2 for strong forces (FN"/k B T > 1). Following the 
discussion about relaxation times for small to moderate 
forces three paragraphs above, we believe that between 
F = and FN" /k B T > 1 the decay for $ R (t) at short 



5 




FIG. 3 (a) Behaviour of <t>(0) — <j>(t) for TV = 100 as a function of t, shown by means of the proxy variable (« (4) (0) - 2 (4) (t)>, 
showing that (f)(0) — <j>(t) reduces to a constant very quickly: F = 1.0 (red) and F = 0.5 (blue). The angular brackets denote an 
average over 672, 000 polymer realizations, (b) Mean time required, for F = 1.0, to unthread a distance s for s = 5, 10, 15, . . . , N: 
N = 100 (red), N = 200 (green), N = 300 (blue). The time-axis corresponding to N = 200 is the true time, for the N = 100 
and N = 300 cases the time axis is divided and multiplied by a factor 2 respectively. This is done in order to show that the 
slope of the curves reduces slowly with increasing N: we obtain, for N = 100 a slope of 0.57, for N = 200 a slope of 0.54, and 
N — 300 a slope of 0.52 at short times. At long times the slope increases for all values of TV: most likely due to the fact that 
the monomer at the pore is too close to the end of the polymer. The solid black line correspond to a slope of 0.52. The angular 
brackets denote an average over 48, 000 polymer realizations. See text for more details. 



times is a combination of these two power laws. Addi- 
tionally, closer inspection of the curves for F — 0.3 and 
F = 0.5 in Fig. |2jb) reveals that the slope is steeper in 
the beginning: perhaps it is an indication that relaxation 
within a blob (corresponding to t l + 2 " ) precedes inter- 
blob rearrangements (corresponding to t^ 1 / 2 ). Neverthe- 
less, at strong forces, the behaviour 



the relation 



Wi(t) ~ t 2 exp(— t/rp), with tf ~ N 2 



(6) 



stands as a witness of the fact that the Flory-like struc- 
ture of the polymer is entirely destroyed. 



IV. THE RELATION BETWEEN <p(t) AND v(t), AND THE 
SCALING BEHAVIOUR OF t n 

A. Relation between the imbalance of chain tension <>(t) and 
the translocation velocity v(t) 

In this subsection we consider the strong force case as it 
is simpler. The moderate to weak force case is discussed 
in Sec. lTVTBI 

So far, we have /iz,(i) ~ i _T +^ exp(— t/rR_ 0Use ) and 
fJ-R(t) ~ t -1 / 2 exp(— £/tf) for strong forces. Since the 
memory effects in the dynamics of the translocating poly- 
mer stem from the power laws, in the absence of symme- 
try between the left and the right segment of the polymer, 
we only need to keep track of the power law of as 
it has a lower exponent than /ii(t). In other words, in 



<t>{t) = <k=o- / dt'\n(t-t')\v(t') 
Jo 



(7) 



we have to use the fact that the power law decay of fi(t) 
behaves ~ i -1 / 2 . Note the absolute value around fJ.(t), 
as the sign of n(t) is negative. 

Equation J7]) can be inverted via Laplace transforma- 
tion, yielding 



v(k) 



(l>t=o _ <Kk) 
fc| M (fc)| | M (fc)| 



(8) 



where k is the Laplace variable representing inverse time. 
Thereafter, using ju(t) - t- 1 / 2 , i.e., n(k) - k- 1 / 2 , and 
Laplace- inverting Eq. , we get 

«(*)= [ dt'(t-t')- 3/2 [<t>t=o-<t>(t)}- (9) 
Jo 



B. Scaling behaviour of tn with N 

In Eq. |(9l), if cj>(t) goes to a constant ^ 4>t=o, then 

v(t)^r^ 2 i.e., s(t)= [ dt' v(t') - t 1 ' 2 , (10) 



where s(t) is the distance unthreaded in time t (|30f). 
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FIG. 4 Collapse of all data in terms of FN V /k B T and 
t n /N 2+ " for F = 0.0, 0.1, 0.2, . . . , 1.0. Inset: the same data 
in log- log plot, the black line corresponds to a slope of —1. 
The rjv values correspond to the median for 1, 024 polymer 
unthreading events. 



In Fig. [3](a) we show the behaviour of [<j>t=o — <j>(t)] 
by means of the proxy variable (z^ 4 ^(0) — z^{t)) for 
strong forces [F = 1.0 (red) and F = 0.5 (blue)], where 
z^\t) is the difference between the Z^(t) values be- 
tween the right and left segment of the polymer, i.e., 



z^\t) = Zf{t) - Z£'(t), Indeed [cj) t=0 - <j>(t)] goes to 



a constant fairly quickly. Following Eq. (flO|) . this yields 
us the scaling s(t) ~ i 1 / 2 for strong forces. The data in 
support of the scaling s(t) ~ t 1 / 2 are shown in Fig. [3jb) , 
for F = 1.0. 

The scaling for the mean unthreading time tn is ob- 
tained from the equation s(tn) — N. For strong forces, 
it is derivable as_7jy ~ N 2 — as shown in Fig. IHa) or 
in earlier works (22 : l23h — from Eq. I|10p if we assume 
that [<j>t=o — 4>{t)] is a constant independent of TV. It 
seems reasonable (and likely!) that a local property like 
[4>t=o — <t>(t)] should be unaffected by the polymer length, 
which is a large-scale property; nevertheless, we have no 
way to argue this theoretically. In the context of using 
Eq. ([Toll to obtain tat ~ iV 2 , it is however useful to note 
that in the scaling sense tn is smaller than (or equal 
to) the time scales in the exponential decay of /XR(t) and 
fJ>L{t), otherwise the power-law behaviour of /z(t) we used 
in Eq. (|1ll0p would not have been applicable for all times 

t < T N . 

The collapse of all the data for the unthreading times 
for several different values of TV and F in terms of the 
variables FN V jk B T and t n /N 2+v , as shown in Fig. H 
indicates that the unthreading time tn can be written in 
a scaling form as 



TN 



f FN V \ 
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(11) 



where g(x) is a scaling function of its argument x. Figure 
H] shows that g(x) ~ 1/x for x ^> 1. Note that 5(0) is 
a constant, in agreement with our earlier result that the 
unthreading time scales as N 2+v in for unbiased trans- 
location |2J), in which case the polymer leaves the pore 
purely due to thermal fluctuations. 

The fact that g(0) is a constant indicates that g(x) has 
to deviate from the 1/x behaviour as x approaches zero. 
From the inset of Fig. [4]we see that g{x) starts to deviate 
from the 1/x behaviour at about x = 4, at which point 
the force is moderate in strength. In Fig. 01 note also 
that tn/N 2 ^ has a higher prefactor close to x = than 
at x = 0. 

A priori, the same analysis (|7lfT0| holds for small to 
moderate forces as well. Nevertheless, whether Eq. ([9]) is 
actually useful in such circumstances is a different mat- 
ter. Indeed, a deeper investigation reveals that <fit=o for 
small to moderate forces can be extremely small. To 
give a feeling for how small <j)t=o can be, we obtained 
(Z^(t = 0)} values for N = 100 for F = 0.0,0.1, ...,1.0 
(not all are plotted in Fig. [p. The (Z {A \t = 0)) val- 
ues for F = 0.1, . . . , 0.3 (approximate x-values 1.5, 3 and 
4.5 in Fig. [4]), corresponding to the right segment of the 
polymer, turned out to be 1.35, 1.36 and 1.38 respectively, 
while (Z^ (t = 0)) corresponding to F = 0.0 (i.e., for the 
left segment of the polymer) turned out to be 1.34. 

Since FN V /ksT is a dimensionless parameter that de- 
scribes the effect of the force on the polymer's dynamics 
compared to the effect of thermal fluctuations, it seems 
logical that if FN V /ksT is slowly reduced, thermal fluc- 
tuations start to dominate over the effect of the force, 
and translocation by the pulling force F starts to resem- 
ble unbiased translocation, i.e., translocation in the ab- 
sence of any external forces ([HI; US H3) • Such a picture is 
manifested by both sides of Eq. © effectively becoming 
zero as suggested in the above paragraph; the equation 
remains valid, but ceases to be useful in practice. 



V. DISCUSSION 

In this paper, we have considered polymer transloca- 
tion pulled through a narrow pore by a force F. We have 
provided a theoretical description of the polymer's dy- 
namics under these conditions, as well as of the scaling 
behaviour of the unthreading time tn for the polymer of 
total length 2N, the time the polymer takes to leave the 
pore. Our theory is supported by high precision com- 
puter simulation data, generated for a three-dimensional 
self-avoiding lattice polymer model. 

At strong forces, i.e., for FN U 3> 1, we have reported 
tn ~ N 2 /F: we have shown that the translocation ve- 
locity v(t) is not constant in time; in fact, the velocity of 
translocation v(t) is shown to behave as £ -1 / 2 , while for 
small to moderate forces the behaviour of v(t) is more 
complicated. At small to moderate forces, satisfying the 
condition FN" < 1, where v « 0.588 is the Flory expo- 
nent for the polymer, we have found that tn is indepen- 
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dent of F , and in agreement with our earlier result for 
unbiased polymer translocation (|24l) scales with polymer 
length as t n ~ N 2+v . 

We have shown that the scaling of v(t) as well as the 
TV-dependent part of tn stem from the dynamics of the 
polymer segments at the immediate vicinity of the pore 
— in particular, the memory effects in the polymer chain 
tension imbalance across the pore. The theoretical anal- 
ysis presented here is based on that of Ref. (24), and 
therefore provides a direct confirmation of the robust- 
ness of the theoretical method presented in Ref. |24T i. 
Additionally, we note that the physical picture provided 
in Refs. |22 : l23h. wherein the scaling arguments for the 
unthreading time involved a constant velocity of trans- 
location (albeit an average one, in light of this work) is 
incomplete. 

It should nevertheless be mentioned that the depen- 
dence of the relevant quantities, such as v(t) or rjy on F is 
beyond the scope of the theoretical description provided 
here. The main reason behind this is that no analyti- 
cal expression has been reported (neither do we have one 
ourselves) for the quantities, such as the polymer chain 
tension, memory kernel etc. for force F. Indeed, the be- 
haviour of the quantities of interest on F is complicated, 
as already manifested by Fig. [TJ and in the absence of a 
theoretical description involving F, numerical investiga- 
tion has remained the only way. Nevertheless, we note 
that the y-axis of Fig . [4] as tn/N 2+1 ' [originating from 
our previous work (24 )]. the z-axis of Fig. |4]as FN" /ksT 
as a measure of the strength of the force in relation to 
thermal fluctuations, and the scaling tn ~ N 2 at strong 
forces automatically imply that tn has to behave ~ 1 /F 
at strong forces. 
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